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ABSTRACT; 

He  present  a  general  approximation  scheme  for  autonomous  FDE's 

with  Lipschitzean  right-hand  side.  Our  approach  is  based  on  approximation 

results  in  semigroup  theory.  The  results  are  applied  to  nonlinear 

autonomous  FDE’s  in  the  state  space  3Rn  x  lA  and  to  linear  autonomous 
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1.  Introduction. 


In  recent  years  one  can  see  considerable  interest  in  approximation  of 
delay  systems  by  ordinary  differential  equations.  Results  of  this  type 
proved  to  be  very  useful  for  the  numerical  solution  of  optimal  control 
problems  and  identification  problems  (see  [2], [5],  [6]).  Avery  successful 
approach  uses  abstract  approximation  results  in  semigroup  theory  and  applies 
these  to  the  semigroups  associated  with  autonomous  delay  systems.  This 
approach  goes  back  to  [4]  (see  also  [5])  where  control  problems  involving 
linear  autonomous  functional -differential  equations  and  approximation  of  the 
state  by  step  functions  are  considered.  A  considerable  improvement  over  the 
scheme  developed  in  [5]  was  obtained  in  [7]  where  for  the  same  class  of  delay 
systems  a  scheme  is  developed  which  allows  approximation  of  the  state  by  spline 
functions.  In  this  paper  we  give  an  abstract  formulation  of  the  scheme 
developed  in  {7]  which  can  be  applied  to  autonomous  delay  systems  with 
globally  Lipschitzean  right-hand  side.  We  also  indicate  the  application  of 
the  linear  version  of  the  scheme  to  autonomous  neutral  delay  systems.  Since 
by  lack  of  space  we  cannot  give  a  detailed  discussion  of  the  relevant  literature, 
we  include  a  rather  complete  list  of  papers  in  the  references  ([1,3,5,7,13,14,15 
16,17,18,19,20,22,23,24])  and  refer  to  the  discussions  given  there  (see  es¬ 
pecially  Section  5  of  [5]). 

2.  An  Approximation  Scheme  for  Semigroups  of  Nonlinear  Transformations. 

Throughout  this  section  X  will  be  a  Hilbert  space  with  inner  product 
(*>*)  and  norm  | • | .  A  family  T(t),  t  0,  of  globally  Lipschitzean 
operators  X  X  is  called  a  semigroup  of  type  to,  to  €  H,  if  the  following 
properties  are  satisfied: 


(ii)  T(t+s)  -  T(t)T(s)  for  all  t,s  >  0. 

(ill)  For  any  x  C  X  the  nap  defined  by  t  •*>  T(t)x  Is 
continuous  on  [0,«°). 

(iv)  |T(t)x  -  T(t)y|  _<  eWt jx-y |  for  all  x,y  €  X  and  t  >  0. 

If  for  a  not  necessarily  single-valued  operator  A  on  X,  the  operator 
(I-XA)-1  is  single-valued.,  and  .defined  on  X  for  X  sufficiently  small,  we 
say  that  T(«)  is  generated  by  A  if  for  any  x  €  X 

T(t)x  “  lim  (1-  —  A)  nx,  t  ^  0, 
rr*=°  n 

uniformly  with  respect  to  t  in  bounded  intervals.  Of  fundamental  importance 
is  the  generation  theorem  by  Crandall-Liggett  ([12]): 

Let  A  be  a  densely  defined  operator  on  X.  If  A  -  wl  is  dissipative 
for  some  to  €  IR  and  if  range  (I-XA)  ■  X  for  sufficiently  small  X, 
then  A  generates  a  semigroup  of  type  u>  on  X. 

If  A  is  single-valued  then  A  -  wl  is  dissipative  if  and  only  if 

a 

(Ax-Ay, x-y)  ^  w|x-y|  for  all  x,y  €  dom  A.  But  note,  that  the  Crandall- 

Liggett  Theorem  is  valid  in  general  Banach  spaces.  Fundamental  for  our  approach 
is  the  following  approximation  theorem  for  semigroups  of  type  u>: 

Theorem  1  ([10]).  Let  A^  N  -  1,2,...,  and  A  be  single-valued  operators 
on  X  such  that  dom  A^,  3  dom  A  for  all  N  and  dom  A  ■  X.  Assume  that  the 


3. 


following  conditions  are  satisfied: 

(i)  There  exists  a  XQ  >  0  such  that 

range  (I-XA)  -  range  (I-XA^)  -  X 
for  N  -  1,2,...  and  all  X  £  (0,X0). 

(ii)  There  exist  real  constants  u^,  N  -  1,2,...,  and  u>  such 
that  the  sequence  {ct^}  is  bounded  above  and  -  u^I 
and  A  -  0)I  are  dissipative. 

.  _(iii)  There  .exists __a.  subset  G  of  dom  A  such  that  (I-XA)G  *  X 

for  X  sufficiently  small  and 


A^x  -*■  Ax  as  N 


for  all  x  €  G. 


Then  A^  and  A  generate  semigroups  T^( •)  and  T(*)  of  type 
and  a),  respectively,  and  for  all  x  €  X 


lim  T„(t)x  »  T(t)x 
***  N 


uniformly  on  bounded  t-intervals. 

This  theorem  follows  Immediately  ffam  Theorems  4.1  and  3.1  in  [10]  and 

* 

is  also  true  in  general  Banach  spaces. 


Definition  1.  Assume  that  A  is  an  operator  on  X.  A  sequence  {X^.P^A^}, 

N  -  1,2,...,  is  called  an  approximation  scheme  related  to  A  if 

(i)  X^  are  finite  dimensional  linear  subspaces  of  dom  A,  N  «  1,2,. 
(ii)  PR:  X  +  Xjj  is  the  orthogonal  projection  onto  X^,  N  ■  1,2,...  . 

(iii)  Aj,  -  PNAPN,  N  -  1,2,...  . 


4. 


Condition  (i)  In  this  definition  is  the  most  restrictive  one.  It  could 
■ean  that  only  the  trivial  approximation  scheme  (0,0,0)  for  all  N  is 
related  to  A.  If  A  generates  a  semigroup  of  some  type  u  it  is  of  interest 
to  know  when  the  operators  A^  generate  semigroups  T^(*)  converging  to  T( •) 

Theorem  2.  Let  A  be  a  densely  defined  single-valued  operator  on  X  such 
that  A  -  wl  is  dissipative  for  some  u  CK  and  range  (I-XA)  =  X  for  X 
sufficiently  small.  Furthermore,  let  {X^,P^ ,A^}  be  an  approximation  scheme 
related  to  A.  Assume  that  the  following  conditions  are  satisfied: 

(1)  11m  P^jx  “  x  for  all  x  6  X. 

(11)  There  exists  a  set  G  c  dom  A  such  that  (I-XA)G  -  X 

for  X  sufficiently  small  and 

11m  AP„x  •  Ax  for  all  xEG, 

W 

(ill)  For  each  N  there  exists  a  X„  >  0  such  that 

N 

range  (I-XA^)  ■  X 
for  all  X  €  (O.X^). 

Then  the  following  is  true: 

a)  Ajj  is  continuous  on  X  and  generates  a  semigroup  of  type  u 
on  X,  N  *  1,2,...  . 

b)  TN(t)XN  c  Xjj  for  tj>0,  N-1,2,...  . 

c)  lim  TM(t)P„x  ■  T(t)x  for  all  x  €  X  uniformly  on  bounded 

M-ko  "  N 
N  ' 

t-lntervals,  where  T(«)  is  the  semigroup  generated  by  A. 
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Proof.  Condition  (11)  of  Theorem  1  is  satisfied  with  *  to.  This  follows 
from 

(AjjX-Ajjy.x-y)  -  (PjjlAPjjX-APjjy]  ,x-y)  = 

2  2 

(APjjX-APjjy.PjjX-Pjjy)  £  wlPNx-PNy|  £u>|x-y| 

* 

for  all  x,y  in  X.  Here  we  have  used  Pjj  “  pjj*  For  x  €  G  we  have 

lV^U  I  PNAPNX_PNA3C  I  +  lPNAx_Axl 
£  |  APNx-Ax  |  +  j  PjjAx-Ax  |  . 

By  conditions  (ii)  and  (i)  both  terms  on  the  right-hand  side  vanish  as  N  -*■  <® 
Thus,  also  condition  (iii)  of  Theorem  1  is  satisfied.  It  only  remains  to 
prove  that  range  (I-XAjj)  =*  X  for  N  ■  1,2,...  and  all  X  £  (0,Xq),  Xq  a 
positive  number  not  dependent  on  N,  and  that  is  continuous  on  X. 

We  have  range  (I-X(AN-u>I))  ■  range  (I-  Ajj)  =  X  for  all  X  with 

0  <  Xd+Xoi)"1  <  XN.  Since  Ajj  -  ul  is  dissipative,  we  immediately  get 
(cf.  [8;  p.73]  or  [9;  p.23])  range  (I-X^-wI))  »  X  for  all  X  >  0,  i.e., 

range  (I-XA^)  •  X  for  X  €  (0,jj).  This  proves  condition  (i)  of  Theorem  1 
with  Xq  -  .  Since  A^  -  col  is  m-disslpative  and  defined  on  all  of  X,  it 

is  also  demicontinuou8,  i.e.,  continuous  from  X  with  the  norm  topology  into 
X  with  the  weak  topology.  But  A^X  c  X^  and  dom  X^  <  00  Imply  that  A^ 

c  Xj,. 

Corollary  1.  If  A  is  linear,  then  the  conclusion  of  Theorem  2  holds  without 
assuming  condition  (iii)  explicitly.  Moreover,  T^(t)  is  given  by 
Tjj(t)  -  exp  Ajjt. 


is  continuous.  Finally,  conclusion  b)  follows  from  A^X^ 


6. 


Proof.  The  conditions  on  A  Imply  (also  in  the  nonlinear  case)  that  A  is 
closed  Ccf .  [8;  p.75]).  Then  also  is  closed  and  defined  on  all  of  X. 

Therefore,  by  the  Closed-Graph  Theorem  AP„  is  a  bounded  linear  operator. 

The  rest  of  the  proof  is  clear. 

In  the  linear  case  we  also  can  replace  the  condition  "range  (I-^A)  *  X 
for  X  sufficiently  small"  by  "A  is  the  infinitesimal  generator  of  a  ^-semi¬ 
group".  This  is  a  consequence  of  the  Lumer-Phillips  Theorem  (see  for  instance 

[21;  p.16]). 

The  following  corollary  covers  a  situation  taylored  for  delay  systems: 

Corollary  2.  Let  A^  be  a  closed  linear  operator  with  dom  Aq  =  X  and  A^ 
be  an  operator  with  dom  A^  3  dom  Aq  such  that 

lAjX-Ajyl  <  L  ||x-y|| 

for  all  x,y  €  dom  Aq,  where  L  >  0  and  ||  •  ||  is  a  norm  on  dom  Aq.  Then 
condition  (Hi)  of  Theorem  2  holds  for  any  approximation  scheme  (x^,PN>  A^ 
related  to  A  »  Ag  +  A^. 

Proof.  The  conditions  on  '■  Aq  Imply  that  A^  q  *  P^AqPjj  is  a  bounded  linear 
operator,  -  | q| .  We  have  to  prove  that  (I-XAg)y  -  x  has  a  solution 

y  for  any  x  6  X  provided  X  is  sufficiently  small.  This  equation  is 
equivalent  to 

y  -  w-UayViv +  -= 

if  X  €  (O.ajj1).  Then  for  all  ylty2  €  x 


7. 


|h(y1)-h<y2)|  <  X(l-aNX)_1|A1PNyi-A1pNy2| 

<  XL(l-a NX)_1  IIP^-P^H  <  OjjXLd^XJ^jy^l 


Here  we  have  used  the  fact  that  two  norms  on  a  finite  dimensional  space  are 
equivalent.  Mote,  that  P^y^  and  PN^2  are  T^e  estimate  shows 

that  h  is  a  contraction  for  X  sufficiently  small. 

Using,  for  instance.  Theorem  II  of  [12]  and  the  fact  that  the  Peano- 
Existence  Theorem  is  valid  in  finite  dimensional  spaces,  we  see  that 
u(t)  *  TN(t)PNx  for  any  x  €  X  is  the  solution  of 


u (t)  -  ^(t) ,  t  ^  0, 
u(0)  -  PNx  , 


which  is  an  initial  value  problem  for  an  ordinary  differential  equation  on 
Xjj.  For  calculations  one  has  to  choose  a  basis  for  and  to  represent 

and  PNx  with  respect  to  this  basis  (see  [7]  and  [15]  for  more  details). 


3.  Applications  to  Functional  Differential  Equations. 
a)  Nonlinear  Equations  in  lRn  x 

We  consider  autonomous  equations 


with  initial  data 


x(t)  -  f(xfc) ,  t  ^  0, 


x(0)  -  n  €Rn,  x(s)  ■  4>(s)  a.e.  on  [-r,0]. 


(3.1) 


(3.2) 
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r  >  0,  where  <j>  €  L2(-r  ,0jIRn) .  For  a  function  x:  [-r,a)  -*-]Rn,  a  >  0,  xfc 

is  defined  by  xt(s)  -  x(t+s),  s  £  [-r,0].  A  solution  of  (3.1),  (3.2)  is  a 

function  x  on  an  interval  [-r,a),  a  >  0,  such  that  (3.2)  holds  and  x(t)  * 

T|  +  f  f(x  )ds  for  t  6  l0,a).  We  assume  that  f  is  a  map  ^^(-r.O^IR®)  -♦■]Rn 
2° 

(Jzf  is  the  linear  space  of  square  integrable  functions)  satisfying 

2  n 

(hi)  For  any  function  x  in  LT(-r ,ci*IR  ),  a  >  0,  which  is  continuous  on 
(Q,a)  t  -*-_f(x£)_  uniquely  defines  a  map  in  L*(0,ajIRn). 

(h2)  There  exist  constants  L  >  0  and  r ^ ,  j  *  0,...,m,  0  «  r^  <  r^  <  ...  <  rm  =>  r, 

such  that  for  any  in  i^2(-r ,0*lRn) 

|f(«  -  fWl  <  U  l  1 4>(— r . )— tp(— r .)  |  +  I  2). 

j=0  J  J  L 

Under  conditions  (hi),  (h2)  problem  (3.1),  (3.2)  has  a  unique  solution 
x(t;n,4>)  on  (-r,00)  for  any  (n,4>)  £  X  =  lRn  *  L2(-r  ,0;IRn) .  The  family  of 
operators  T(-)  defined  by  T(t)(n,<J>)  “  (x(t;n,<j>)  ,xt(n,<l>)) ,  t  ^  0,  (n,<l>)  £  X, 
is  a  strongly  continuous  semigroup  of  globally  Lipschitzean  operators  on  X 
(cf.  (17]).  X  is  a  Hilbert  space  with  norm  |  (n ,4)  |  =  ( | n | 2  +  [  |<J>(s)  |2g(s)ds)^2 

and  corresponding  inner  product,  where  |*|  denotes  the  Euclidean  norm  in  ]R 
and  the  weighting  function  g  is  defined  by 

g(s)  -  m-j+1  for  s  £  (-r^.-r^j),  j  -  l,...,m. 

Define  the  operator  A  by 

doe  A  -  {  (4(0)  ,$)|4>  €  W1»2(-r,0jB*1)}, 

A(*<0),$)  -  (f(*)^),  +  €  W1,2(-r ,0jRn) . 


i.'m j  ^  .  nui il  u1. "i  lwjmm 


Theorem  3.  Let  {X^,P^,A^}  be  an  approximation  scheme  related  to  A  such 
that 

(i)  lim  P  x  =  x  for  all  x  €  X. 

N-ko  N 

(ii)  There  exists  an  integer  k  1  such  that  £or4>€  Ck(-r,OjIRn) 

$N(0)  -  <f>(0)  and  ^  $  in  L2(-r,0jIRn) 

as  N  Here  <J>N  £  WX,2(-r,0;lRn)  is  defined  by 

--  -  pN(<Kt)) ,<)))  »  (<t>N(0),(j)N)  £  dom  A. 

Then,  for  each  N,  Ajj  generates  a  semigroup  T^(*)  on  X  such  that  for  all 
x  €  X 

lim  T  (t)PNx  =  T(t)x 

N-k» 

uniformly  on  bounded  t-intervals. 

Indication  of  proof  -(a  complete  representation  of  the  results 
concerning  equation  (3.1)  will  appear  elsewhere  [15]). 

Step  1;  We  have  dom  A  ■  X.  For  the  set  G  in  condition  (ii)  of  Theorem  2 
we  take  G  =>  (  (<j>(0)  ,<j>)  |  <\>  £  Ck(-r  ,0  jRn)  } .  Using  (h2)  one  shows  that  there  exists 
an  u  6  3R  such  that  A  -  ul  is  dissipative  in  X.  The  idea  to  introduce  a 
weighting  function  is  due  to  G.F.  Webb  [25]  and  is  essential  here.  In  order 
to  verify  the  range  condition  on  A  one  observes  that  (I-XA)  (<j>(0)  ,<j>)  * 

(n,<!0  €  X  is  equivalent  to  $(s)  *  e8^(0)  "  ^  |  e(s-T)  A^(T)dT f  s  €  [-r,0], 

and  4>(0)  -  h(<J>(0)),  where  h(a)  -  n  +  Xf(e‘/Xa  -  £  j  e^^VOdT)  for 
a  £]Rn.  Using  (h2)  one  shows  that  h  is  a  contraction  on  K.  for  X 


sufficiently  small. 
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Step  2 1  By  definition  of  A  we  see  that  AP^x  -*•  Ax  in  D  is  equivalent  to 

f  (4>jj)  "*  f  C<J>)  and  4>N  -*■  4>  in  L^(-r,0jIRn) .  But  f(<J>jj)  "*■  f  (4>)  follows  from 

(ii)  and  (h2) .  Therefore,  condition  (il)  of  Theorem  2  is  satisfied. 

Step  3:  In  order  to  apply  Corollary  2  we  write  A  -  Aq  +  A^,  where 

AqOKOM)  =  (0,*)  and  A1(<K0),<|>)  -  <f  (<J>)  ,0)  for  $  €  W1»2(-r,0'»n) .  Then 

A^  is  a  densely  defined  and  closed  linear  operator  (see  [17]).  (h2)  implies 

that-  At  satisfies  the  condition  in  Corollary  2  with  ||  (4> CO)  ,<J>>  ||  =  sup  |<Ks)| 

-r<s<0 

Thus,  condition  (iii)  of  Theorem  2  is  satisfied. 


Step  4 :  The  semigroup  generated  by  A  is  indeed  T(«)  (cf.  [17]). 


The  conclusion  of  Theorem  3  in  the  case  of  spline  approximation  was  also 
obtained  by  H.T.  Banks  in  [3]  by  a  different  method  under  the  additional 
assumption  that  f  is  differentiable. 


b)  Equations  of  Neutral  Type  in  IRn  x  L^(-r ,QdRn) 

Let  D  and  L  be  bounded  linear  functionals  C(—r,0jRn)  IRn 


defined  by 


D0J>) 

L(4>) 


m  rO 

<p(0)  -  l  B.<|>(-r.)  -  B(s)<j>(s)ds, 

j-1  3  3  '  -r 

m  f0 

l  A.<j»(-r.)  +  A(s)$(s)ds , 
j-0  3  3  J-r 


where  0  ■  rn  <  r,  <  . . ,  <  r  ■  r.  A.,B,  are  n  *  n  matrices  and  A(*)>B(’) 
ul  m  j  j 

are  in  L^(-r,0j]RnXn) . 


As  before  we  take  X  **  3Rn  x  L2(-r,0;IRn) .  Following  an  idea  given  in 
[11]  we  define  the  operator  A  by 

don  A*  (  (D(<J>)  ,<j>)  |  €  W^,2(-r,OpRn)}, 

A(D(4>)  ,<t>)  =  (L(<J>)  ,4>) ,  4>  €  W1,2(-r,0^Rn)  . 


X  will  be  supplied  with  the  norm  |(n,<J>)|  .  where  the  weighting  function  g 

8 

now  is  also  dependent  on  D: 


g(s) 


m-j+l+ct^ ( s+r j _^)  for  s  €  (-r^  .-r^) ,  3  *  1, 


m(nH-2)  ■  <2  .  . 

<  “  -  I®4i  »  J  “  1»  •  •  •  »m* 

3  "  rj  rj-l  J 


Theorem  A. 

a)  A  generates  a  ^-semigroup  T<  * )  of  bounded  linear  operators  on  X. 

b)  Let  {Xjj,Pjj,Ajj}  be  an  approximation  scheme  related  to  A  satisfying: 

(i)  lim  P  x  *  x  for  all  x  6  X. 

N-*» 

(ii)  There  exists  an  integer  k  _>  1  such  that  for 

<f>  €  Ck(-r,0-IRn)  <f>N(0)  -*■  4»(0)  and  <frN  -*•  4>  in  L2(-r,0-IRn) 

as  N  -*•  where  ^  €  W*,2(-r ,0jlR.n)  is  defined  by 

pn(d(4>),$)  -  (d(*n),4»n). 


Then,  for  each  N,  A^  is  a  bounded  linear  operator  which  generates  the  (^-semi¬ 
group  T^(t)  ■  exp  Ajjt,  t  ^  0,  and  for  all  x  €  X 

lim  T^OPjjX  -  T(t)x 
IH® 


uniformly  on  bounded  t -intervals. 
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The  proof  of  this  theorem  proceeds  in  similar  steps  as  that  for  Theorem 
3  and  can  be  found  In  [14].  If  z  ■  (D($),$)  with  $  6  C(-r,0dRn)  then 
T(t)z  *  (D(x^)  ,x^),  t  _>  0,  where  x(t)  is  the  solution  of  the  neutral  equation 

ft  D(xfc)  -  L(xfc),  t  >_  0, 

with  x(s)  -  <Ks),  s  €  [-r,0].  Theorem  4  provides  an  approximation  of 
t  -*■  (D(xt),xt),  t  >_  0,  in  X.  If  we  could  prove  that  TN(t)PNz  -*■  T(t)z  also 
with  respect  to  the  sup-norm,  then  we  would  get  an  approximation  of  x(t) 
itself  (see  [14],  Section  6  for  details  and  preliminary  results). 

1  O  1% 

c)  Equations  of  Neutral  Type  in  Vr*  (-r,03R  ) 

Let  the  operators  D  and  L  be  asunder  b)  and  put  D(<J>)  -  <J>(0)  -  D($) . 
The  Cauchy  problem 

•  m  r0  ^ 

x(t)  -  L(xt)  +  l  B  x(t-r,)  +  B(s)x(t+s)ds  -  L(xfc)  +  D(xt>,  t  >_  0, 
j-1  3  3  '-r  ** 

x(s)  -  4>(s),  s  €  [-r,0j, 

has  for  any  <J>  €  W^,2(-r,0 jRn)  a  unique  solution  x(t ;4>)  on  [-r,00).  We  take 
X  *  W*,2(-r,0;IRn)  and  define  the  family  T(*)  of  operators  by 

T(t)<|>  -  xt (<p),  t  >  0,  $  c 

T ( • )  is  a  Co-8emigroup  of  bounded  linear  operators  on  X  with  infinitesimal 
generator  A  given  by 


. . 


dom  A  -  £  W2,2(-r,0 jIRn)  and  4(0)  =  D(4>  +  L(<j>)}, 

A4  *  4>  <p  6  dom  A. 

Note,  that  dom  A  depends  on  D  and  L.  We  supply  X  with  the  norm 

■  |<f>L  ■  (|<J>(0)|2+  [  |^<s)  j2g(s)ds)1/2,  where  g  is  as  under  b)  with 

8  J-r 

a*  <  -  ~  2m . -  |B,|2,  j  =  1 . m. 

3  ~  3 

Then  there  exists  an  w  €  K  such  that  A  -  ojI  is  dissipative  in  X  (cf. 

1161). 

Theorem  5.  Let  {X^.P^A^}  be  an  approximation  scheme  related  to  A  such 
that : 

CD  lim  P  4  «  4  for  all  4  €  X. 

N-**> 

(ii)  There  exists  an  integer  k  >_  2  such  that  for  4  E  C^C-r.OpR11) 
4n<0)  -*■  4(0)  and  4N  +  4  ki  L  (-r,0;tttn)  as  N  ■*■  «>, 
where  4N  €  W2,2(-r,0jRn)  is  defined  by  PH4  -  4N< 

Then  each  A^  is  a  bounded  linear  operator  which  generates  the  (^-semigroup 
TN(t)  -  exp  A^t  on  X  and  for  all  4  €  X 

lim  TN(t)PN4  -  T(t)4 

N-*eo 

uniformly  on  bounded  t-intervals. 

For  a  proof  of  this  theorem  see  [16].  Theorem  5  provides,  for  initial 
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data  in  W1,2(-r,0*IRn) ,  an  approximation  of  xfc  with  respect  to  the  sup-norm 
•  2 

and  of  x£  with  respect  to  the  L  -norm.  On  the  other  hand,  the  subspaces 
Xjj  and  projection  depend  on  the  right-hand  side  of  the  equation  even 

if  D(<|>)  -  4>(0). 

4.  Spline  Approximation 

In  this  section  we  indicate  that  subspaces  of  spline  functions  can  be 
used  in  order  to  get  concrete  algorithms  using  the  schemes  presented  in  Section 
3. 

In  the  situation  of  Section  3,  a)  and  b)  we  can  take  any  sequence  of 


N 

subdivisions  of  [-r,0]  with  mesh  points  t  ^ ,  j  «*  1,...,^,  such  that 
max|tj-tj_^|/min|tj-tj_^|  £  (5  <  °°.  In  case  of  Section  3,  a)  we  define  Xjj  to 

be  the  subspace  of  all  elements  (4(0),<j>)  such  that  <f>  is  a  first  order, 

N 

cubic  or  cubic  Hermite  spline  with  knots  at  the  t ^ ,  respectively.  He  have 
dim  X^  •  (k^+l)n,  (k^+SJn  or  2(kjj+l)n,  respectively.  Similarly  in  the  case 
of  Section  3,  b)  is  the  space  of  all  elements  (D(<j>)  ,<J>)  with  <|>  as 

above. 

In  the  situation  of  Section  3,  c)  we  have  to  satisfy  the  boundary 

condition  <K0)  -  D(<f>)  +  L(<J>)  and  to  guarantee  more  smoothness.  One  has  to 

take  subdivisions  of  [-r,0]  such  that  the  points  -r j ,  j  ■  0,...,m,  are 

meshpoints  and  restrict  D  and  L.  In  the  case  of  cubic  splines  we  have  to 

assume  »  ...  -  Bm_^  •  0,  B(s)  =  A(s)  =  0.  Then  is  the  subspace  of 

all  cubic  splines  with  knots  at  the  meshpoints  t,  satisfying 

.  m  .  ^ 

♦(0)  ■  l  A.<K-r,)  +  B  (-r) .  In  the  case  of  cubic  Hermite  splines  we  assume 
jiO  J  J  m 


B(s)  =  A(s)  =  0  and  take  to  be  the  subspace  of  all  cubic  Hermite 

N  ®  tn 

splines  with  knots  at  the  t”  such  that  <K0)  =  £  A.4>(-r.)  +  I  B  <K-r  ) . 

J  j»0  J  3  j«l  J  J 

For  details,  explicit  calculation  of  bases  for  X^,  matrix  representa- 
tions  of  A^  and  numerical  examples,  see  [7],  [14]  and  [16]. 
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